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FOREWORD 

The  Defense  Communications  Engineering  Center  (DCEC)  Technical 
Notes  (TN's)  are  published  to  inform  interested  members  of  the  defense 
community  regarding  technical  activities  of  the  Center,  completed  and  in 
progress.  They  are  intended  to  stimulate  thinking  and  encourage  infor- 
mation exchange;  but  they  do  not  represent  an  approved  position  of  policy 
of  DCEC,  and  should  not  be  used  as  authoritative  guidance  for  related 
planning  and/or  further  action. 

Comments  or  technical  inquiries  concerning  this  document  are  welcome, 
and  should  be  directed  to: 

Director 

Defense  Communications  Engineering  Center 

1860  Wiehle  Avenue 

Res ton,  Virginia  22090 
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I.  INTRODUCTION 


In  transmitting  blocks  of  data  through  a communications  system,  ensuring 
that  the  system  is  in  synchronization  is  vital.  Each  frame  or*  block  of 
data  contains  a header  of  information.  These  headers  are  checked  to  ensure 
that  the  system  is  still  in  sync.  There  are  two  possible  ways  that  the 
header  may  not  be  correct:  first,  some  bits  may  be  transmitted  wrong  or 
secondly,  the  system  could  be  out  of  sync.  Once  sync  has  been  lost,  some 
corrective  action  has  to  be  taken, during  which  the  system  may  be  down.  In 
order  to  ensure  that  the  system  is  not  taken  down  for  bursts  of  transmission 
errors  an  up/down  counter  is  used  to  monitor  its  behavior. 

The  counter  has  N+1  slots  labelled  0,1,..., N and  a pointer  that  is  moved 
up  or  down  according  to  whether  the  header  is  correct  or  not.  Usually,  if 
the  header  is  correct  the  pointer  is  moved  up  one  slot  - unless  it  is  at  N in 
which  case  it  remains  there.  The  number  of  slots  that  the  pointer  is  moved 
down  is  the  subject  of  this  investigation.  We  analyze  the  case  where  the 
pointer  is  moved  down  K(£N)  slots  when  the  header  is  wrong.  We  assume  that  if 
the  pointer  is  within  K slots  of  zero,  then  the  pointer  is  moved  to  zero 
when  an  incorrect  header  appears.  Once  at  zero,  the  pointer  remains  there 
and  the  system  is  assumed  to  be  out  of  sync.  The  focus  of  the  study  is  the 
random  variable  which  describes  the  elapsed  time  from  starting  at  slot  N until 
the  pointer  is  at  zero  for  the  first  time.  This  random  variable  represents 
the  time  until  checking  the  system  to  see  if  it  is  still  in  sync  once  the 
counter  has  been  reset. 

If  p is  the  probability  an  arriving  header  is  wrong,  and  q=l-p  that  it 
is  right;  then  the  system  can  be  modelled  as  a Markov  chain  and  the 
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analysis  can  be  carried  out  in  a straightforward  manner.  In  section  II, 
this  is  done;  section  III  contains  some  numerical  examples.  Finally, 
section  IV  gives  some  concluding  statements. 

This  problem  was  generated  because  of  conversations  the  author  has  had 
with  Dr.  Stan  Lorens  of  R710  and  Mr.  Walter  Roehr  of  R740.  Dr.  Lorens  pointed 
out  that  the  Naval  Research  Lab  (NRL)  has  written  a specification  for  an  Ad- 
vanced Narrowband  Digital  Voice  Terminal  (ANDVT)  [1],  The  linear  predictive 
coder  in  this  terminal  will  be  capable  of  extracting  frame  sync.  The  frame 
synchronization  shall  consist  of  two  modes;  sync  acquisition  and  sync  main- 
tenance. The  sync  maintenance  function  is  the  one  that  is  analyzed  in  this 
technical  note.  For  the  ANDVT  it  consists  of  a counter  numbered  0 to  15, 

When  a correct  sync  bit  is  received,  the  counter  is  incremented  by  one.  If 
the  check  bit  is  incorrect  it  is  decremented  by  two.  Once  the  counter 
reaches  zero,  sync  has  to  be  reacquired.  So  the  ANDVT  sync  maintenance  is  a 
special  case  {N=15,  K=2)  of  our  general  counter. 

Dr,  Lorens  stated  that  to  the  best  of  his  knowledge  no  analysis  had  been 
done  to  determine  if  the  NRL  specification  was  satisfactory,  Mr.  Roehr  added 
that  it  would  be  desirable  to  have  the  capability  to  work  at  more  general  parameter 
than  just  N=15  and  K=2  as  in  the  NRL  specification.  These  conversations  have 
resulted  in  the  mathematical  analysis  and  computer  implementation  of  these 
results  presented  in  this  technical  note.  The  author  is  grateful  to  Dr.  Lorens 
and  Mr.  Roehr  for  bringing  this  problem  to  his  attention. 

As  a closing  cormient  to  this  section  we  point  a related  piece  of  work  by 
Smith  [2].  He  considered  the  more  general  problem  of  acquiring  and  maintaining 
sync.  In  this  paper  a more  extensive  reference  list  of  related  work  on  this 
problem  is  given. 
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II.  MATHEMATICAL  ANALYSIS 


As  we  mentioned  in  the  preceding  section,  the  system  can  be  modelled  as  a 
Markov  chain.  For  a discussion  of  these  types  of  stochastic  processes,  see 
[3],  [4],  [5]  or  [6].  We  will  use  some  definitions  from  Markov  chains. 

The  reader  should  refer  to  the  texts  cited  above  for  an  explanation  of  un- 
familiar terms. 

Let  be  the  random  variable  representing  the  position  of  the  pointer  just 
after  the  header  was  observed  for  the  nth  time.  We  adopt  the  convention 
that  the  observation  of  the  header  and  the  subsequent  changing  of  the  placement 
of  the  pointer  occur  simultaneously.  The  random  variables  X^,  n=l,2,..., 
can  take  on  the  values  0,1,..., N.  Since  the  probability  that  the  header  is 
correct,  q,  does  not  depend  on  n, the  value  X^  takes  on  only  depends  on  X^_i 
and  as  such  forms  a Markov  chain.  In  Markov  chains  the  most  important  element 
which  governs  how  the  chains  progress  in  time  is  called  the  probability  transi- 
tion matrix,  usually  denoted  by  P.  The  i,jth  entry  of  this  matrix  gives  the 
probability  that  X^=j  given  that  X^_i=i  for  all  i and  j.  For  the  Markov  chain 
being  considered  here  the  probability  transition  matrix,  P,  is  given  by 


N-1  0 

N 0 


P 0 

p 0 


0 


0 q 
0 q 


The  zero  state  is  called  an  absorbing  state  in  that  once  this  Markov 
chain  goes  there.it  never  leaves.  If  q=l.  the  state  N would  also  be  absorbing, 
but  in  the  application  to  up/down  counters,  q is  probably  close  to  but  not 
equal  to  1 . 

Now  we  are  in  position  to  define  the  random  variable  that  we  wish  to 

study.  Suppose  at  some  instant,  denoted  by  0,  we  observe  the  system  and  the 

pointer  is  pointing  at  the  ith  slot,  i.e.,  XQ=i . We  want  to  study  the  be- 
havior of  the  system  as  it  goes  to  zero  for  the  first  time.  We  define  Y^-  to  be 
the  random  variable  representing  the  time  (in  number  of  observations  of  the 
header)  until  the  pointer  is  at  zero  for  the  first  time,  given  that  the  initial 

observation  of  the  pointer  is  at  the  ith  slot.  We  are  really  interested  in 

Y^^  but  need  the  more  general  definition  for  the  development.  Note  that  Y,^ 
represents  the  time  until  checking  for  sync  once  the  counter  has  been  reset. 

Let  us  define 


fi^n= 


. . . , j'O  I XQ=i 


We  are  interested  in  studying  fj^  p’ probability  of  the  pointer  being 
at  zero  for  the  first  time  at  the  nth  observation  of  the  header  given  that  it 
was  initially  set  at  the  Nth  slot.  For  a general  K,p,  and  N the  exact  solution 


for  f.,  is  difficult;  some  special  cases  will  be  discussed  in  section  IV,  but 
N,n 

a simple  recursion  relation  allows  one  to  iteratively  solve  for  fj^^  For 
n=l  we  have  (fg  g=l,  f^  q=0  i>l) 


^ H 0 here)  we  have 

i+1, n-1  ■*'P^O,n-l 

i=l,2,... 

i+l,n-l'^P^i-K,n-l 

i=K+l,K+2 

N.n-1'^P^N-K,n-1 

i = N . 

(4) 


The  development  of  equation  (4)  is  as  follows:  if  the  initial  observation  of 
the  pointer  is  at  the  ith  slot  and  we  want  to  go  to  zero  at  the  nth  observa- 
tion for  the  first  time.we  consider  what  happens  at  the  first  observation. 

For  i=l,2,...,K  we  can  go  to  the  (i+l)st  slot  with  probability  q (header 
correct);  once  at  the  (i+l)st  step  we  have  to  go  to  zero  for  the  first  time  in 
n-1  observations.  The  latter  quantity  is  f.^.|  by  the  homogeneity  of  the 
Markov  chain.  With  probability  p we  would  go  to  zero  from  these  slots  but 
n^2  in  equation  (4),  and  since  we  stay  at  zero  there  is  no  contribution  to 

f.  „ for  this  event, 

1 ,n 

For  i=K+l ,K+2, . . . ,N-1  we  can  move  up  one  slot  with  probability  q to  i+1 
or  drop  back  K to  i-K  with  probability  p.and  then  go  to  zero  in  (n-1)  observa- 
tions. Finally,  for  i=N  we  stay  at  slot  N with  probability  q or  drop  back 
to  N-K  with  probability  p.  The  form  of  equation  (4)  makes  it  ideally  suitable 
for  implementation  on  the  computer. 

The  exact  solution  fj^  ^ for  a general  K is  difficult,  but  a form  of  solution 
be  obtained  via  generating  functions.  Define  for  i=l,2,...,N  and  |z|;<l  (iIq(z)=1) 


ni(z) 


00 


(5) 


5 


then  from  the  first  part  of  equation  (4)  we  have 

n^(z)=qzn.^^(z)+pzi  (6) 

the  last  term  on  the  right  hand  side  of  equation  (6)  comes  from  f . . For 

* » I 

i=K+l,...,N-l 

n.(z)=qzn.^^(z)+pzn._^(z)  (7) 

and  finally  for  i=N  in  equation  (4),  one  gets 

(l-qz)nj^^(z)=pzn,,^_^(z).  (8) 

The  solution  technique  for  n^(z)  is  a standard  one  for  these  types  of 
equations  and  follows  the  development  given  in  [5].  Equation  (7)  suggests  a 
trial  solution  of  the  form  (z)=x’ (z) . Substituting  this  form  of  solution  in 
equation  (7)  one  gets 

qzx'^'^^(z)-A'^(z)+pz=0.  (9) 

Let  X^(z),  r=l ,2, . . . ,K+1 , be  the  roots  forthis  equation,  then  the  general  solu- 
tion to  equation  (7)  is  given  by 

K+1  . 

n.(z)=  z a (z)x  (z),  (10) 

1 r=l 

where  a^(z),  r=l ,2, , . . ,K+1 , are  to  be  determined. 

These  K+1  unknowns  may  be  found  from  the  K equations  given  in  equation  (6) 
and  the  one  given  by  equation  (8).  Using  these  equations  and  equation  (10)  one 
gets  the  following  K+1  equations  in  a^(z);  for  i=l,2,...,K 

I a (z)[xj.(z)(l-qzx  (z))]=pz  (11) 

r=l  ^ 
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and 


2:  o‘-(z){(l-qz)xJ|(z)-p2xJ|"'^(z)}  = 0.  (12) 

r=i 

Thus,  equations  (9),  (11)  and  (12),  can  be  used  to  solve  the  problem, 
ror  practical  application,  this  solution  becomes  cumbersome  although  the  re- 
sults for  K=1  provides  some  insight  and  is  now  given.  For  K=1,  there  are 
two  roots  of  equation  (9);  they  are 


and 


X (2)  . 

' o 


2qz 


,(2)  . 

2qz 


(13) 


(14) 


Since 

n^.(z)=a^(z)Ai(z)+a2(z)>^2^z), 

using  equation  (11)  for  K=1,  one  finds 


1 = (2)+a2(z) . (15) 

Using  this  result  in  equation  (12),  a2(z)  is  given  by 

(l-qz)x']'(z)-p2X,^"\z) 

ct2(z)  = ! (^°1 

pz(x^^-''(z)-x!J"''(z))+(l-qz)(x^(z)-x^^(z)) 

These  results  may  be  used  to  find  an  analytic  expression  for  the  expected  value 
of  Y^,  E(Y^),  That  is, 

E(Yf^)=n'(l)=  alP^A-^Np^q-p)  ^ 

^ ^ p^2p-l)(q-p) 

For  a general  K,  a simple  expression  for  even  E(y|^)  is  difficult  to  obtain. 
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Although  one  may  use  equations  (6),  (7),  and  (8)  to  recursively  find  all 
the  moments  of  we  find  the  mean,  E(Yfj),  and  variance,  Var(Y|^),  of  Yj^ 

from  n,^(z).  Using  the  basic  definitions  we  have 
E(Y^)  = n’(l) 


V3r(Y|^)  = n|^(l  )+ii|^(l  )-Ti^(1 ) . 

Let  n .=iT^  (1 ),  then  using  equations  (6)  and  (7)  we  have 

E.  = 1+qE.^^  ; i=l,2,...,K  08) 

and 

E^  = l+qE.^^+pE._,^  : i=K+l  ,K+2, . . . ,N-1 . (19) 

Note  we  have  used  the  fact  that  n.(l)=l.  From  equations  (18)  and  (19)  one 
sees  that  E^  for  i=2,3,...,N  can  be  expressed  in  terms  of  E-j  so  let  us  define 

E.  = C(i,l)E^+C(i,2);  (20) 

that  is,  E-  is  a constant, C(i ,1 ),  times  E^  plus  another  constant  C(i,2), 

with  C(l,l)  and  C(l,2)=0.  Using  equation  (18)  we  have 


rC(i-l,l) 

qC(i,l)=| 

Lc(i-l,l)-pC(i-l-K, 


i=2,3,... ,K+1 


i=K+2,.., ,N 


C(i-1,2)-1 


rni-i, 

0.2)  = 

tc(i-i. 


2)-l-pC(i-l-K,2) 


i=2,3,...,K+l 


i=K+2,...,N 


From  equation  (8)  we  have 


pE^  = 1+pEm  ,,, 
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which  can  be  used  with  equations  (19)  and  (2Q)  to  find 


^ ^ 1+P(C(N-K.2)-C(N.2))  ^ (24) 

^ P(C(NJ)-C(N-K.l))  ’ 

Thus,  one  can  recursively  evaluate  C(i,l)  and  C(i,2),  using  equation  (24)  to 
find  and  then  set 

s E(Y^)  = C(N,1)E^+C(N,2). 

II 

In  order  to  find  Var(Y|^)  one  has  to  find  n|^(l);  this  can  be  done  in  the 

I 

same  manner  we  found  ilj^(l).  For  i=l,2,...,N,  let 
= Il!,'(l)  = D(i,l)F^+D(i,2), 

obviously,  D(i,l)=l  and  D(i,2)=0.  Proceeding  as  before  we  have 


rD(i-l.l) 

qD(i,l)  = 

LD(i-l,l)-pD(i-l-K,l) 


:i=2,3, , . . ,K+1 
:i=K+2 N 


and 


qD(i,2) 


D(i-l,2)-2qE. 

D(i-l,2)-pD(i-l-K,2)-Z^ 


:i=2,3,...,K+l 

:i=K+2,...,N 


where  Z^.=2qE^+2pE^._.|_l^.  Again,  using  equation  (8)  we  have 


(25) 


(26) 


^ p[D(N-K,2)-D(N,2)]  + 2qE,^+2pEfj.l^ 

^ p[D(N,l)-D(N-K,l)] 

and  F^i  = D(N,1  )F^+D(N,2) . 


(27) 


1^ 
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We  have  given  a recursive  method  for  easily  computing  E(Y|^),  Var(Yj^)  and 
the  complete  probability  density  function  of  Y|yj.  In  the  next  section  we 
give  some  numerical  examples  using  these  results. 


III.  SOME  NUMERICAL  EXAMPLES 


In  this  section  we  give  three  numerical  examples  using  the  results  pre- 
sented in  section  II.  The  first  is  given  in  Table  I;  the  probability  density  func- 
tion of  the  number  of  observations  of  the  header  until  the  pointer  goes  from 
N=15  to  zero  for  the  first  time  is  given.  For  that  table,  K=2  (i.e.,  the 
NRL  specification)  and  p was  set  equal  to  .6,  .7,  .8,  and  .9.  The 
reason  for  the  large  values  of  p is  to  ensure  that  non-zero  values  ^ 

are  basically  given  when  n^40.  Since  K=2  we  have  f^^  ^=0  for  n<_7;  that  is, 
it  takes  at  least  8 observations  for  the  system  to  drop  to  zero.  It  is  inter- 
esting to  point  out  the  oscillatory  nature  of  f^g  ^ for  a fixed  value  of  p. 
Consider  p=.8,  then  f^g  g=. 26844,  f^g  ^q=. 05369  and  f^g  ^^=.16106;  thus 
these  probabilities  depend  on  K and  how  the  sample  function  of  Y-jg  may  proceed 
to  zero.  ' 

I 

Tables  II  and  III  give  values  of  the  mean  and  variance  of  the  time  until  ’ 

checking  for  out  sync  for  N=15  and  fixed  different  values  of  p and  K.  One  of 
the  most  interesting  facts  that  can  be  seen  in  these  figures  is  that  ' 

E(Yl5)^:^Var(Y^g),  (28)  | 

(where  means  approximately)for  p small.  This  fact  was  not  true  for  the 
parameters  used  in  Table  I.  The  fact  that  the  parameters  selected  for 
Tables  II  and  III  are  more  representative  of  real  world  values  tends  to  make 

j 

the  result  given  by  equation  (28)  more  promising  for  future  work.  This  idea  j 
will  be  discussed  in  the  next  section.  I 
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Another  interesting  and  important  fact  that  can  be  seen  is  that  the  I 

performance  of  this  up/down  counter  is  extremely  sensitive  to  the  values  of  | 

) 

p and  K.  From  Table  II,  for  K=2  when  p changes  from  .01  to  .5,  the  expected  I 

1 5 ^ 

time  to  check  to  see  if  the  system  is  out  of  sync  changes  from  1.09x10  to  \ 

275.  If  a header  is  arriving  every  .0225  seconds,  then  this  number  corres-  | 

ponds  to  checking  for  sync  either  every  18.7x10®  years  or  every  6.1875  seconds. 
Similarly,  the  p fixed  and  K varying  we  have  the  same  sort  of  results;  see 
p=.05  in  Table  II.  These  results  indicate  that  setting  the  specification  at 
a given  value  of  K is  extremely  critical  and  should  be  studied  in  more  depth. 
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TABLE  I.  PROBABILITY  DENSITY  FUNCTION  FOR  TIME  UNTIL  CHECKING  IF 
OUT  OF  SYNC  (N=I5,  K=2) 


^N.n 

n 

P 

.6 

.7 

.8 

.9 

8 

.01680 

.16777 

.43047 

9 

.05375 

. 26844 

.34437 

10 

.02150 

.04151 

.05369 

.03444 

n 

.05375 

.11414 

. 1 61 06 

.11192 

12 

.09374 

.15627 

.15247 

.05458 

13 

.03956 

.04950 

.03221 

.00577 

14 

.06525 

.08748 

.06073 

.01156 

15 

.09392 

.10130 

.04967 

.00496 

16 

.04176 

. 03384 

.01107 

.00055 

17 

.05748 

.04909 

.01689 

.00089 

18 

.07520 

.05252 

.01296 

.00036 

19 

.03495 

.01835 

.00302 

.00004 

20 

.04348 

.02378 

. 00407 

.00006 

21 

.05373 

.02432 

.00302 

.00002 

22 

.02592 

. 00883 

.00073 

.00000 

23 

.03023 

.01003 

.00091 

24 

.03596 

.01056 

.00066 

25 

.01790 

.00396 

.00017 

26 

.01998 

.00453 

.00019 

27 

.02312 

.00441 

.00014 

28 

.01183 

.00170 

.00004 

29 

.01278 

.00187 

. 00004 

30 

.01449 

.00130 

.00003 

31 

.00759 

.00071 

. 00001 

32 

.00800 

.00076 

.00001 

33 

. 00892 

.00072 

.00001 

34 

.00477 

.00029 

.00000 

1 

35 

.00493 

.00030 

/ 

36 

.00343 

.00028 

37 

.00296 

.00012 

38 

.00301 

.00012 

39 

.00328 

.00011 

\ 

/ 

N 

/ 

40 

.00182 

.00005 

.00000 

.00000 

m 

-< 

18.209 

13.603 

10.880 

9.109 

Var(Y^5) 

52.539 

18.828 

7.080 

2.219 

13 
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TABLE  II.  EXPECTED  TIME  UNTIL  CHECKING  IF  OUT 
OF  SYNC  FOR  N=15* 


.01 

to 

o 

0.1 

0.2 

0.3 

Mtm 

0.5  ■ 

■ 

8.87x10^^ 

1.78x10^^ 

2.9x10^^ 

2.39x10^ 

1.45x10® 

6.48x10° 

2.40x10^ 

2 

1.09x10^^ 

1.70x10^ 

3.27x10® 

6.40x10^ 

2.36x10^ 

5.47x10^ 

2.75x10^ 

3 

6.62x10^ 

7.68x10^ 

1.17x10^ 

2.28x10^ 

4.69x10^ 

2.25x10^ 

1.46x10^ 

D 

2.07x10^ 

1.97x10^ 

9.08x10^ 

6.81x10^ 

2.46x10^ 

1.43x10^ 

1.01x10^ 

5 

6.94x10^ 

3 

2.64x10-^ 

2.46x10^ 

3.7x10^ 

1.68x10^ 

1.07x10^ 

7.8x10° 

6 

9.67x10"^ 

7.62x10^ 

1.14x10^ 

2.51x10^ 

1.28x10^ 

8.54x10° 

6.4x10° 

■ 

3.74x10^ 

3.89x10^ 

7.35x10^ 

2.01x10^ 

1.11x10^ 

7.8x10° 

6.07x10° 

8 

4.68x10^ 

1.71x10^ 

4.60x10^ 

1.51x10^ 

n 

8.6x10'^ 

5.95x10° 

4.51x10° 

9 

2.55x10^ 

1.16x10^ 

3.57x10^ 

1.28x10^ 

7.55x10° 

5.33x10° 

4.13x10° 

10 

1.78x10^ 

9.22x10^ 

3.04x10^ 

1.16x10^ 

7.08x10° 

5.11x10° 

11 

1.39x10^ 

7.84x10^ 

2.73x10^ 

1.09x10^ 

6.86x10° 

5.04x10° 

4.01x10° 

12 

1.14x10^ 

6.95x10^ 

2.53x10^ 

1.06x10^ 

6.76x10° 

5.02x10° 

4.00x10° 

13 

9.79x10^ 

6.34x10^ 

2.39x10^ 

1.04x10^ 

6.71x10° 

5.01x10° 

4.00x10° 

14 

8.61x102 

5.9x10^ 

2.30x10^ 

1.02x10^ 

6.69x10° 

5.00x10° 

4.00x10° 

15 

1.00x10^ 

2.0x10^ 

1.0x10^ 

5.00x10^ 

3.33x10° 



2.50x10° 

o 

o 

X 

o 

o 

CNJ 

I *This  time  is  the  number  of  sync  maintenance  checks  that  have  to  be  made.  To  con- 

I vert  this  to  time  one  has  to  multiply  this  number  by  the  expected  time  to 

I i transmit  a frame. 

I 1 

.1 


TABLE  III.  VARIANCE  OF  TIME  UNTIL  CHECKING  IF  OUT  OF 
SYNC  FOR  N=15 


1^ 

.01 

.05 

0.1 

0.2 

0.3 

0.4 

0.5 

1 

59 

7.86x10^^ 

3.17x10®® 

8.5x10^® 

5.69x10^® 

2.09x10^® 

4.14x10^ 

3.83x10^ 

2 

1.21x10^® 

2.87x10^® 

1.4x10^® 

4.07x10^ 

4.57x10^ 

1.34x10® 

1.91x10® 

3 

4.39x10^^ 

5.9x10^^ 

1.4x10® 

4.57x10® 

1.2x10® 

1.69x10® 

4.79x10^ 

B 

4.30x10^^ 

3.87x10® 

8.1x10® 

3.39x10® 

2.81x10® 

6.56x10^ 

2.38x10^ 

5 

4.8x10^^ 

6.91x10® 

5.65x10"^ 

8.94x10® 

1.24x10® 

3.66x10^ 

1 51x10^ 

6 

9,35x10^ 

5.71x10® 

1.14x10^ 

3.92x10® 

7.3x10^ 

2.42x10^ 

1.02x10^ 

■ 

1.39x10^ 

1.46x10® 

4,52x10® 

2.27x10® 

4.77x10^ 

1.65x10^ 

7.2x10® 

8 

2.19x10^ 

2.8x10^ 

1.78x10® 

1.44x10® 

3.79x10^ 

1.53x10^ 

7.4x10® 

1 

9 

6.51x10® 

1.28x10^ 

1.05x10® 

1.04x10® 

2.91x10^ 

1.16x10^ 

5.4x10® 

10 

3.16x10® 

7.85x10® 

7.39x10^ 

8.17x10^ 

2.34x10^ 

9.4x10® 

4.5x10® 

11 

1.91x10® 

5.54x10® 

5.71x10® 

6.79x10^ 

2.01x10^ 

8.4x10® 

4.15x10® 

12 

1.29x10® 

4.25x10® 

4.69x10® 

5.89x10^ 

1.81x10^ 

7.9x10® 

4.05x10® 

13 

9.47x10® 

3.45x10® 

4.01x10® 

5.29x10^ 

1.70x10^ 

7.7x10® 

4.01x10® 

14 

7.31x10® 

2.91x10® 

3.53x10® 

4.89x10^ 

1.63x10^ 

7.6x10® 

4.00x10® 

15 

9.9x10^ 

3.80x10® 

9.0x10® 

2.0x10^ 

7.77x10^ 

3.8x10® 

2.0x10® 

15 


IV.  CONaUSIONS 


[ 


We  have  given  a mathematical  analysis  for  the  operation  of  an  up/down 
counter  that  is  used  to  monitor  the  sync  behavior  of  data  in  a telecommunica- 
tions system.  We  close  this  technical  note  with  a few  remarks  concerning  the 
work  and  some  possible  extensions. 

First,  from  a historical  point  of  view  the  formulation  of  the  problem  pre- 
sented here  is  closely  related  to  the  classical  gambler's  ruin  problem.  This 
problem  is  thoroughly  discussed  in  Feller  [3].  In  that  problem  K=1  and  the 
state  N is  also  an  absorbing  state.  We  have  not  searched  the  literature  on 
the  gambler's  ruin  problem  thoroughly  enough  to  ensure  that  more  work  on  our 
problem  has  not  been  done.  If  further  work  in  this  area  is  to  be  done  one  should 
look  at  the  literature  more  closely. 

The  results  that  we  gave  were  in  the  form  of  recursive  relationships.  As 
such  they  are  easily  implemented  on  the  computer.  Two  problems  arise  with  this 
type  of  solution:  first,  since  no  analytic  expression  (except  for  K=1 ) was 
given  for  fj^  E(Yn^),  or  Var(Y|^),  the  exact  dependency  on  N,  K,  and  p can  only 
be  ascertained  by  drawing  conclusions  from  the  numerical  examples.  It  would 
be  desirable  to  have  an  analytic  expression  so  that  more  concrete  conclusions 
can  be  drawn;  although  this  appears  to  be  a formidable  task. 

The  second  problem  with  the  recursive  solution  is  one  of  numerical 

accuracy.  The  smallest  value  of  p that  we  have  considered  in  the  numerical 

examples  is  p=  0.01,  for  K<^5  and  N=15;  this  gave  values  of  ^(Yig)  and  Var(Yig) 

7 14 

on  the  order  of  10  and  10  or  greater.  Nimibers  of  this  magnitude  are  subject 
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to  the  numerical  accuracies  of  the  computer  one  is  dealing  with.  If  closed 
form  analytic  expressions  were  given, this  problem  might  be  avoided  as  was  done 
for  the  K=1  case. 

For  N=K,  a little  thought  will  convince  the  reader  that 

f^^^  = :n=l,2 (29) 

from  which 

E(Yfj)=  P'^  (30) 

and 

Var(Y^)=  qp-^  (31) 

These  results  along  with  equation  (17)  allow  us  to  give  some  bounds  on 
E(Yj^)  as  a function  of  K.  For  N and  p fixed  E(Yf^)  is  minimum  when  K=N  and 
maximum  when  K=l.  Thus,  equation  (17)  is  an  upper  bound  and  equation  (30)  a 
lower  bound  on  the  expected  number  of  observations  until  the  pointer  is  at 
zero  for  the  first  time. 

We  remarked  in  section  III  that  for  the  values  of  p used  in  Tables  II 
and  III,  the  mean  time  until  out  of  sync  was  approximately  equal  to  the 
square  root  of  the  variance.  Since  some  numerical  problems  could  be  en- 
countered in  evaluating  the  performance  of  the  counter  for  p's  smaller  than 
.001,  some  sort  of  approximation  for  these  cases  should  be  considered.  One 
possible  approach  is  along  the  lines  of  a diffusion  approximation  as  pre- 
sented by  Feller  [2]  for  the  gambler's  ruin  problem.  It  would  be  interesting 
to  see  if  that  development  could  be  extended  to  our  problem. 
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Finally,  it  should  be  reemphasized  that  the  table  given  in  section  III 
clearly  shows  that  the  setting  of  K for  a fixed  value  of  p is  extremely 
critical.  Just  changing  the  value  of  K one  slot  can  make  several  orders  of 
magnitude  difference  in  the  results  one  gets.  Extreme  care  should  be  taken 
when  fixing  K.  For  instance,  for  the  ANVDT  we  had  N=15,  K=2.  If  the 

probability  of  the  check  bit  being  received  wrong  (p)  changes  from  .01  to  .5 

this  corresponds  to  having  to  reacquire  sync  every  18.7x10^  years  to  every 
6.1875  seconds. 

It  is  hoped  that  the  results  presented  in  this  paper  may  be  used  to 
provide  some  insight  into  the  performance  of  up/down  counters  used  to  ensure 
the  system  is  maintaining  sync. 


t 
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